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Lipid membranes regulate the flow of materials and information 
between cells and their organelles. Further, lipid composition and 
morphology can play a key role in regulating a variety of biological 
processes. For example, viral uptake, plasma membrane tension 
regulation, and the formation of caveolae all require the creation 
and control of groups of lipids that adopt speciflc morphologies. 
In this paper, we use a simplifled model mixture of lipids and 
cholesterol to examine the interplay between lipid phase-separation 
and bilayer morphology. We observe and theoretically analyze 
three main features of phase-separated giant unilamellar vesicles. 
First, by tracking the motion of 'dimpled' domains, we measure 
repulsive, elastic interactions that create short-range translational 
and orientational order, leading to a stable distribution of domain 
sizes, and hence maintaining lateral heterogeneity on relatively 
short length scales and long time scales. Second, we examine 
the transition to 'budded' domain morphologies, showing that 
the transition is size-selective, and has two kinetic regimes, as 
revealed by a calculated phase diagram. Finally, using observations 
of the interactions between dimpled and budded domains, wc 
build a theoretical framework with an elastic model that maps 
the free energies and allowed transitions in domain morphology 
upon coalescence, to serve as an interpretive tool for understanding 
the algebra of domain morphology. In all three cases, the two 
major factors that regulate domain morphology and morphological 
transitions are the domain size and membrane tension. 

Cellular membranes are a complex mixture of lipids, mem- 
brane proteins, and small molecules {e.g. sterols) [1] [2]. 
The membrane serves mainly as a chemical barrier and sub- 
strate for membrane proteins that are responsible for regu- 
lating the movement of materials and information across the 
membrane. However, there are a host of important tasks 
that require a change in membrane morphology, such as endo- 
and exo-cytosis |3] |4] , vesicular trafficking from the endoplas- 
mic reticulum and Golgi apparatus to the plasma membrane 
[5], and the regulation of tension in the plasma membrane 
[6]. While the role of proteins cannot be ignored in these in- 
stances {e.g. clathrin, COPI, COPII, caveolin, SNAREs, actin) 
[Zl[ll|9l[ini[IIl[l2l[I3[ll[l5l[l6],the lipid composition and 
bilayer morphology of the membrane play an important part 
[5l[T7l[l8l[lll[20l[lTl[22]. With that in mind, our goals in this 
paper are to examine how lipids in a model multi-component 
membrane spatially organize, how this organization relates to 
membrane morphology, or specifically membrane mechanics, 
and to examine how transitions in membrane morphology are 
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regulated by bilayer mechanical properties. 

In vitro studies have conclusively shown that lipids are 
capable of lateral self-organization |23l 1241 125) . facilitated 
by the structure of their hydrophobic regions and the pres- 
ence of intercalated sterols. Saturated lipids and cholesterol 
are sequestered from the membrane mix to form 'lipid rafts' 
that serve as platforms for signaling and material transport 
across the membrane, with sizes ranging from ~ 50 — 500 nm 
[mSlElEZlIlllIIlESlEIlElIMl ED- in addition to their 
unique chemical properties and protein-specific interactions, 
lipid rafts are mechanical entities in a thermal environment, 
and as such, our analysis focuses on continuum and statistical 
mechanics of phase-separated bilayers. 

The remainder of the paper is organized as follows. The first 
section examines how dimpled domains, which exert repulsive 
forces on each other, spatially organize themselves into ordered 
structures as their areal density increases. Specifically, we 
measure the potential of mean force and orientational order, 
finding that domains exhibit orientational and translational 
order on length scales much larger than the domains them- 
selves. The second section studies the transition to a spherical 
'budded' morphology [351 136) . Using a mechanical model that 
combines effects from bending, line tension and membrane 
tension, we predict and observe size-selective budding transi- 
tions on the surface of giant unilamellar vesicles, and derive a 
phase diagram for the budding transition. The final section 
considers the three lipid domain morphologies - flat, dimpled 
and budded - and constructs a set of transition rules that 
dictate the resultant morphology resulting from coalescence of 
two domains. 

Spatial Organization of Dimpled Domains 

Our analysis begins by viewing the bilayer as a mechanical en- 
tity, endowed with a resistance to bending [37], quantified by 
a bending modulus {nb) with units of energy; a resistance to 
stretch under applied membrane tension (r) [37], with units 
of energy per unit area; and in the case where more than one 
lipid phase is present, an energetic cost at the phase boundary, 
quantified by an energy per unit length (7) [38U39| . For a given 
domain size, the line tension between the two phases competes 
with the applied membrane tension and bending stiffness to 
yield morphologies that reduce the overall elastic free energy. 
In particular, the bending stiffness and membrane tension both 
favor a flat domain morphology. Conversely, the line tension 
prefers any morphology (in three dimensions) that reduces the 
phase boundary length. A natural length-scale, over which per- 
turbations in the membrane disappear, is established by the 
bending stiffness and membrane tension, given by A = \J K^jr. 
Comparing this 'elastic decay length' with domain size indi- 
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cates the set of possible domain morphologies. If the domain 
size is on the order of, or smaller than, the elastic decay length, 
flat and dimpled morphologies arise; while domains larger than 
A generally give rise to the budded morphology |35l 1361 1401 141j . 
This rule of thumb is based on the fact that bilayer deforma- 
tions from dimpling are concentrated within a few elastic decay 
lengths of the phase boundary, while domain budding energet- 
ics are governed by basis shapes much larger than A. 

Dimpled domains are characterized by a dome-like shape 
with finite slope at the boundary between the two lipid phases, 
as shown in Figs. [T]and[2jb). In previous work we found a 
distinct flat-to-dimpled transition, regulated by line tension 
and domain area [42) . where if all the elastic parameters are 
constant in time, the domain is either flat or dimpled, but 
cannot make transitions between those states {i.e. there is no 
coexistence regime). Hence, domains below a critical size lie 
flat, and if all other membrane properties are constant, the 
only way domains transition from flat to dimpled is by coa- 
lescing to form larger domains that are above a critical size. 
An important outcome of domain dimpling is the emergence of 
a membrane-mediated repulsive interaction between domains 
that tends to inhibit coalescence. Intuitively, the origin of this 
force is that dimples deform their surrounding membrane, but 
this deformation decays back to an unperturbed state within a 
few elastic decay lengths of the phase boundary. Two domains 
that are within a few elastic decay lengths of each other have 
overlapping deformed regions, and thus the total elastic free 
energy depends on the distance between the domains, leading 
to a net repulsion. A relatively simple mechanical model and 
previous measurements show that the interaction between two 
dimpled domains can be approximated by the pair potential 
V{r) oc e^""/^, where r is the center-to-center separation be- 
tween the domains |42) . 

This repulsive interaction arrests coalescence, and hence sig- 
nificantly affects the evolution of domain sizes in a phase- 
separated membrane. For a simple physical model, in the ab- 
sence of any interaction, domains would diffuse [33] and coa- 
lesce at a rate such that domain size scales as t^^^ [44l |45] |46] . 
This sets the time-scale for full phase separation on a typical 
giant unilamellar vesicle (GUV) at ~ 1 minute (see Fig. [U^c) 
and [23]). However, viewing repulsive interactions as an ener- 
getic barrier to domain growth, and given the measured barrier 
height of ~ SfcsT (with fcs = 1.38 x 10"^^ J/K and T ~ 300 K) 
[42j , the rate of domain coalescence slows by the Arrhenius fac- 
tor e~"' — 0.007. A clear example of the difference in the rate of 
domain growth, with and without elastic interactions, is shown 
in Fig.[2jd and c), respectively, and [47]. Thus elastic repulsion 
is a plausible mechanism by which lipid lateral heterogeneity 
could be maintained on the hour-long time scale required for 
a cell to recycle (and hence partially homogenize) the plasma 
membrane [IH]. Alternative schemes have been proposed that 
balance continuous rates of membrane recycling and domain 
coalescence to yield a stable domain size distribution |36) . 

We examined the role these elastic interactions play in the 
spatial organization of lipid domains. Given that all the do- 
mains mutually repel each other, as the areal density of do- 
mains increases, the arrangement of domains that minimizes 
the elastic free energy takes on distinctly hexagonal order, so 
as to maximize the separation between all domains. Indeed, the 



arrangement of repulsive bodies on a sphere is a well studied 
problem in physics |49l 1501 151) , and a dominant feature of such 
systems is the emergence of hexagonal and translational order. 
We measured the strength of this organizing effect by track- 
ing the thermal motion of dimples on the surface of GUVs, 
and calculating the radial distribution function (see 'Materi- 
als and Methods'). For time-courses that have no coalescence 
events, the vesicle and its domains are in quasistatic equilib- 
rium, thus the negative natural logarithm of the radial distri- 
bution function is a measure of the potential of mean force 
between domains. Our previous theoretical and experimental 
work showed that the elastic interaction between domains at 
low areal density is well approximated by a pair potential of 
the form V{r) oc e~^^^ [42] . As the domain areal density in- 
creases, the domains adopt spatial orientations that maximize 
their mean spacing. This can be understood in terms of the 
free energy of the entire group of domains. If a domain deviates 
from this spatial arrangement, the sum of the elastic interac- 
tion energy with its neighboring domains increases, providing 
a mild restoring force to its original position. Thus the po- 
tential of mean force develops energy wells, up to ~ 2fc_Br in 
depth, that confine domains to a well-defined spacing, as shown 
in Fig. [3Ib-f). It should be noted that such a restoring force 
can arise from the combined pair repulsion of the hexagonally- 
arranged neighboring domains, and does not necessarily mean 
that there are attractive, non-pairwise interactions. 

The mutual repulsion and resulting energetic confinement 
lead to an effective lattice constant that depends on domain 
size and packing density. For example, domains may exhibit a 
well-defined spacing for first, second (Fig.[3jb)), third (Fig.[3jc- 
d)), and fourth (Fig. [3{e-f)) nearest-neighbors, corresponding 
to a correlation in the position of domains over a few microns. 
Interestingly, this means that by forming dimpled domains, the 
motion of individual lipids can be correlated on length scales 
up to ~ 10'' times larger than the size of an individual lipid. As 
indicated by the exponential decay of these 'ringing' potentials 
(see Fig. [3fb-f)), the length scale over which this correlation in 
motion exists is limited by both the relatively low strength of 
pair repulsion (relative to ftsT) and the dispersion of domain 
sizes. 

In the picture that emerges, lipid domains exhibit a tran- 
sition similar to condensation in a liquid-gas system. At the 
lowest areal densities, the motion of domains is analogous to 
a 'gas' of particles that occasionally have repulsive pairwise 
interactions. As the domain areal density increases a 'con- 
densed' phase of domains emerges, identified by its transla- 
tional and orientational order. As the areal density of dimpled 
domains increases, the system exhibits three qualitative effects: 
i) the lattice constant decreases, as it must, to accommodate 
more domains per unit area; ii) the effective confinement grows 
stronger because the membrane in between the domains is 
more severely deformed by the closer packing; iii) hexagonal 
order clearly emerges, as shown by the characteristic peaks in 
the time-averaged Fourier transforms of domain positions in 
Fig- EJb-f) (inset). The time-averaged Fourier transform is the 
arithmetic mean of the Fourier transforms of domain positions 
from each image in a data set, with the peaks corresponding to 
hexagonal order somewhat 'smeared' by the rotational diffusion 
of the entire group of domains. 
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To quantify the lattice constant and correlation length of in- 
teracting dimpled domains, we added a phenomenological cor- 
rection term to the previously mentioned pair potential to ac- 
count for interactions between multiple domains, such that the 
total potential of mean force has the form 

Vr,,{r) = a^e-^ + a2e-^ Jo{2-K{r - r,)/\z) , (1) 

where Jo is the O^'-order Bessel Function of the first kind, and 
ai, Ai, a2, A2, To and A3 are fit parameters. While this equa- 
tion offers little insight into the underlying physics of densely 
interacting domains, the description does an excellent job of 
capturing the observed features of interactions between mul- 
tiple dimpled domains, as demonstrated in Fig. [31 Using this 
formula, we extracted the correlation length (A2) and lattice 
constant (A3), whose ratio, A2/A3, is a measure of the trans- 
lational order in the system, which is shown to increase with 
domain density. As domain areal density increases, the elas- 
tic free energy confines domains to adopt a well-defined mean 
spacing with hexagonal order. Thus the motion of domains 
is correlated over multiple layers of neighboring domains, as 
shown in Fig. |4l 

The dimpled domains that exhibit this behavior arise in sit- 
uations where the tension is low and the elastic decay length 
is longer than the domain size. In the regime where the elastic 
decay length is short compared to domain size, 'budded' do- 
mains emerge as a morphology with distinct transition rules 
and interactions. 

The Budding Transition 

Similar to the analysis of domain dimpling; bending, mem- 
brane tension, line tension, and domain size all play a role in 
the transition to a budded domain morphology. Many ener- 
getic models have been proposed that describe morphological 
changes which result in budding and other more complex mor- 
phologies go] ED [53 El EH [56]. One of the simplest models, 
and yet most reconcilable with experiment, is the 'spherical' 
budding model. This model has its foundations in classical 
'sessile' droplet wetting theory [57] . and posits that the domain 
is, at all times, a section of a sphere [6] [35]. We will recapit- 
ulate this model here, and explore some of its implications for 
our experiments. This model ignores deformations near the 
phase boundary, and cannot capture the existence of the dim- 
pled state, but is a reasonable model to employ in the regime 
where the elastic decay length is smaller than the domain size. 

The budding domain is characterized by a wrapping angle 
6, where 6 = corresponds to a flat domain and 6 = tt corre- 
sponds to the encapsulation of a small volume by a spherical 
bud, as shown in Fig. [TJc) and [SJa) . The bending energy of 
a budding domain is calculated as a fraction of the bending 
energy of a sphere, given by 

Gbond = 2nt J (H^ Cof AA (2) 

where H — 1/R is the mean curvature, R is the radius of 
curvature of the domain, Sttki, is the bending energy of a sphere, 
A is the domain area, and Co is the spontaneous curvature of the 
domain, which, for simplicity, we assume is zero. A constant 



energy C, that does not depend on domain shape, is omitted. 
As the domain becomes more spherical, the areal footprint of 
the domain shrinks, as shown in Fig. [5{a), and work must be 
done against the applied membrane tension, given by 

Gtan. ^ -TTT{Rsm6f. (3) 

The driving force for budding is the reduction of phase- 
boundary line tension, provided by 

Giino = 727ri?sin6'. (4) 

Finally, for all reasonable membrane tensions, the domain area 
is conserved during any change in morphology, and is given by 

^ = 27r_R^(l - cose*). (5) 

This constraint equation links A and R, allowing us to elimi- 
nate R from the total free energy, G = Gbcnd + Gtona + Gii„o. 
After some rearrangement, the total free energy can be written 
in a compact form, 

^ t r- / 1 + COS e I+COS6I 

G = 47rKb|xv"Y — g:;;^ — g:;;^ — + i - cosg j, (6) 

^ ■ V ■ ' ^ ^-^^ — ^ bonding 

where the dimensionless area, a = A/X^, and dimensionless 
line tension, x = 7A/Kt,, emerge as the regulators of domain 
budding. The stable morphologies are those found at energy 
minima, given by dG/dO — 0, with only flat {9 = 0) or budded 
{9 — tt) morphologies satisfying this equation (in the absence 
of spontaneous curvature) . Figure [5]^b) shows the free energy 
of budding as a function of wrapping angle 6 and the line ten- 
sion X- From this plot, one can readily see that there are two 
special values of the line tension; the first, shown in blue, is 
where the free energy difference between the flat and budded 
states equals zero, but an energy barrier exists between them. 
The second special value of line tension, shown in red, is where 
the energy barrier between flat and budded morphologies dis- 
appears, and budding becomes a spontaneous process. This 
graphical analysis primes us to calculate the budding phase di- 
agram. From the solutions for the energy minima, it can be 
shown that the phase diagram has three regions, as shown in 
Fig. [S] i) for certain values of a and x both flat and bud- 
ded domains are stable (coexistence), but flat domains have a 
lower elastic free energy; ii) in an adjacent region, both mor- 
phologies are stable (coexistence), but buds have a lower elastic 
free energy; iii) in the remaining region only buds are stable 
(single-phase) . The boundary between the regions of the phase 
diagram that have two stable morphologies (coexistence) versus 
one stable morphology (single-phase) is given by the inflection 
point {d'^G/d9^)\^_^^ = 0, which defines the line tension 

„ [tt fa ,„, 
Xbud = 8J- + J-, (7) 
\j a \ n 

above which only buds are stable, or alternatively stated, there 
is no energy barrier to the budding process (see Fig. [6}. Given 
that X is a constant material parameter for constant tension, 
this equation specifies a size range over which spontaneous do- 
main budding will occur, 

J(x-yZ^)'<a< J(x+yF^)'. (8) 
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Thus budding, and in particular spontaneous budding, is a 
size-selective process that can only occur if x > 4\/2. Mem- 
brane tension and line tension can be estimated by measuring 
size-selective spontaneous budding on the surface of a phase- 
separated vesicle. In a few instances, we were able to cap- 
ture the onset of size-selective budding, though as a function 
of initial conditions and timing, this proves particularly dif- 
ficult. Sample temperature is a coarse knob that allows us 
to change the state of tension on the vesicle surface. Though 
the exact value of the thermal area expansion coefficient for 
bilayers varies with composition, a good approximate value 
in the temperature range of interest is Ccxp — 5 x 10"^-?^"^ 
[58II59II50] . In Fig.[7{a-c), the temperature is increased slightly 
(from ~ 18C to ~ 20 C, see 'Materials and Methods'), in- 
creasing vesicle area by approximately 1% while maintaining 
the enclosed volume, thus lowering the tension and driving the 
system into the spontaneous budding regime. The average size 
of budding domains is r = 0.93 ± 0.18 ^m. Using eqn. [8l and 
taking ftt = 25fesT as a nominal value for the bending mod- 
ulus of a domain [241 1371 1611 162) , we can solve the equations 
defined by the upper and lower bound to find the line ten- 
sion and membrane tension. From this analysis, we estimate 
r ~ 2.4 X 10~* kBT/nrn^ and 7 ~ 0.45 fcijT/nm. Using the 
tension and our assumption of bending modulus, we can also 
calculate the elastic decay length and dimensionless domain 
size to find A ~ 320 nm and a ~ 26. This membrane tension, 
which is within the range set by typical free vesicle [241 I37[ 
and unstressed plasma membrane experiments [631 164] . sets 
the dimensionless domain area larger than one, and hence sug- 
gests that the spherical budding model is a good approxima- 
tion. This estimate of line tension is consistent with previous 
measurements [24II65| . and quantitatively matches results from 
our previous work [42| . 

As a prelude to the calculation of allowed morphological 
transitions, we note that for the morphology of a domain to 
move from one region of the phase diagram to another, the 
domain must either change size via coalescence, or there must 
be a change in membrane tension which affects both a and x- 
On the phase diagram in Fig. |6l horizontal lines would corre- 
spond to increasing domain area, and the dashed trajectories 
are increasing membrane tension with fixed domain area. The 
key fact is that, except within a region very near the phase 
boundary, the free energy difference between the flat and the 
fully budded states is much larger than knT, as is the energy 
barrier between those states {e.g. Fig. [5j. Thus, from an equi- 
librium statistical mechanics perspective, a budding domain 
can be approximated as a two state system, with the sponta- 
neous budding regime included within the budded state. Thus, 
it makes sense to impose the thermodynamic requirement that 
the free energy difference between morphological states be neg- 
ative for a transition to be allowed, i.e. G\e=TT — G\e^o < if 
going from flat to budded. This amounts to describing budding 
with a two-state model where 

AG/_6 = G\e=^ - G|e=o = ttrs (p^ - 2xp + 8) , (9) 

and p = \/a/n is the dimensionless domain radius. Figures 
(TJ^a-c) and[8ljb-d) show the two states of budding on the surface 
of phase-separated vesicles. If we consider p, a measure of do- 
main size, as an independent variable, then the single control- 



parameter X dictates whether the thermodynamic condition 
AGf^b < has been met for a particular domain size. If the 
dimensionless line tension is below the critical value Xc = 2\/2, 
defined by AG/^i, = 0, the budding transition is forbidden for 
all domain sizes. If x > Xc, budding is allowed, thought not 
necessarily spontaneous, within the size range given by 

P = X±Vx^-8, (10) 

as demonstrated in Fig. [7jf)- This size range always includes 
the range specified by eqn. (8] because spontaneous budding 
always has a negative free energy. 

The 'Algebra' of Morphology 

With an understanding of the conditions under which a domain 
transitions from flat to dimpled [42], and dimpled to budded, 
we are in position to calculate the change in free energy when 
domains of different morphologies coalesce. On a short enough 
time-scale, coalescence only occurs between two domains at a 
time, and hence we can think of the coarsening behavior of 
a phase-separated membrane as many such binary coalescence 
events happening in succession. The purpose of this section is 
to begin to build a framework for understanding how domain 
morphology and coalescence work in concert to affect the mor- 
phological evolution of a phase-separated membrane. In partic- 
ular, we calculate the allowed, resultant morphology when two 
domains, each of a distinct morphology, coalesce. The change 
in free energy associated with a change in domain morphology, 
from flat to dimpled, or dimpled to budded, is much greater 
than ksT, and hence, like the budding analysis of the previous 
section, for a particular transition to be allowed, we demand 
that the change in free energy be negative. Furthermore, the 
large reduction in line energy upon coalescence (compared to 
ksT) means that, in general, coalescence is irreversible, and 
hence after each coalescence event the system is presumed to be 
in a unique quasistatic equilibrium state, with a unique mem- 
brane tension. The use of these transitions rules must then be 
considered in the context of these unique states, that is, tran- 
sitions involving domains of a particular size that were allowed 
before a coalescence event might be prohibited afterward, or 
vica versa. 

Let us denote transitions that involve flat domains with the 
letter /, dimpled domains with the letter d, and budded do- 
mains with letter b, such that, for instance, a flat domain coa- 
lescing with a dimpled domain to yield a budded domain would 
be denoted by fd —> b. There are six possible binary coales- 
cence events: //, fd, fb, dd, db, and bb; each resulting in a sin- 
gle domain of either /, d or 6 morphology. Thus at the onset, 
there are a total of 18 possible morphological transitions, how- 
ever, not all of them are thermodynamically allowed. Specif- 
ically, any domain whose size is greater than the critical size 
required for dimpling cannot adopt a flat morphology as there 
is no flat-dimple coexistence, hence only the // f transi- 
tion can end with an / domain (see Fig.[5fa)). This eliminates 
flve of the six possibilities that end with an / domain. The 
only other transition that can be eliminated immediately is 
// b, because the coalescence of two flat domains must first 
go through the dimpled state. 

This leaves twelve possible morphological transitions, as 
shown in Fig.[9ja-^). For simplicity we will assume the domains 
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have no spontaneous curvature (though this is straightforward 
to incorporate |35II66) ). The free energy change associated with 
each of these twelve transitions is calculated by knowing the 
free energy change associated with three simpler transitions, 
namely the / — » d, / — > 6, and ff^f transitions. Of these, 
the f ^ b transition free energy was discussed in the previous 
section, and the f ^ d transition free energy is a complicated 
function discussed at length in [32], though we note the impor- 
tant fact that AGf^d{ai + 02) < AG/^d(Q:i) + AG/^d(Q2) 
if both domains are above the critical size for dimpling, or in 
words, the free energy of domain dimpling as a function of do- 
main area grows faster than linearly. The scheme we are about 
to build is a valid frame work for understanding energy based 
transitions because we know that changes in free energy, when 
moving along a reaction coordinate, are additive. 

The transition of two flat domains coalescing to yield another 
flat domain is the most fundamental transition, as shown in 
Fig. [21 and can be calculated as the difference in the line tension 
energy between the initial and final states given by 

AG//-^/(ai,a2) = 2^/7rKi,x [V^i -f 02 — s/a^ — , (11) 

where ai and 02 are the dimensionless areas of the two domains 
and AG f f-,f{ai, a^) < —ksT for all domain areas of one lipid 
or more. This situation, depicted by Fig. [H^a) and shown ex- 
perimentally in Fig. [2jc), is encountered at high membrane 
tension, when domains are too small to dimple before and af- 
ter coalescence. 

Using the three basic transitions, we now address the re- 
maining eleven transitions in detail. The next transition we 
consider is two flat domains, each too small to dimple on their 
own, coalescing to form a domain large enough to dimple, as 
depicted in Fig. [9jb). The transition free energy is given by 

AG//_,d(ai, 02) = AG//_,/(ai, 02) + AG/^d(Qi + 02), 

(12) 

and is negative as long as oi + 0L2 is greater than the critical 
area required for dimpling. 

The next transition is a flat and dimpled domain coalesc- 
ing to form a dimpled domain, as depicted in Fig. [U^c). The 
transition free energy is given by 

AG/d-.d(ai,02) = AG//^d(ai,02) - AG/_d(a2). (13) 

No definitive statement about the resultant morphology after 
coalescence of a fiat and dimpled domain can be made, be- 
cause the free energy of this transition must be compared to 
the closely related transition of a flat and dimpled domain co- 
alescing to form a budded domain, to determine which has a 
greater reduction in free energy. This related transition, de- 
picted in Fig. [Hl^h), has the transition free energy 

AG/d-.6(ai,Q!2) = AG//_/(ai,a2) - AG/-,d(a2) (14) 

+AG/^6(ai -f a2). 

Which of these two transitions, fd —> d or fd b, domi- 
nates depends on which has a greater reduction in free en- 
ergy. Comparing eqns. [131 and 1141 asking which transition has 
the greater reduction in free energy is simply asking whether 
AG/^d(Qi-|-a2) > AG/^i,(ai -|-a2) or vice versa. This energy 
balance between the f ^ d and f b transitions determines 



the outcome of all of the subsequent binary transitions as well, 
though we will rigorously show this for the remaining cases. 
Because this energetic comparison crops up so often, we will 
simply refer to it as the 'bud-dimple energy balance.' 

The next transition is two dimpled domains coalescing to 
yield a dimpled domain, as depicted in Fig. [Hd) and shown in 
Fig. [HJa) . The transition free energy is given by 

AGdd^d(Qi, 02) = AGfd^d{cei,a2) - AG/_d(a2). (15) 

Again, we must consider a closely related transition, namely the 
coalescence of two dimpled domains yielding a budded domain, 
as depicted in Fig. [5{i) , with transition free energy 

AGdd^b{cei,a2) = AGdd^d(Qi, 02) - AG/^d(Qi + 02) (16) 

-|-AG/^6(ai + Q2). 

Comparing these two related transitions, dd d and dd — > b, 
we see that the dominant transition is determined by the bud- 
dimple energy balance. 

The next transition is a flat and a budded domain coalesc- 
ing to form a dimpled domain, as depicted in Fig. [9|g). The 
transition free energy is given by 

AG/6-.d (01,02) = AG//_d(ai,a2) - AG/_6(a2). (17) 

The related transition, where a flat and budded domain coa- 
lesce to form a budded domain, as depicted in Fig. [9j£), has 
the transition free energy 

AG/6^6(ai, 02) = AG/6^d(ai, 02) - AG/-.d(ai + Q2) (18) 

+AGf^b{cil + C12)- 

Comparing these two related transitions, fb—>d and fb b, 
we see that the dominant transition is determined by the bud- 
dimple energy balance. 

The next transition is a budded and a dimpled domain coa- 
lescing to form a budded domain, as depicted in Fig. [9jj) and 
shown in Fig. [8jc-e) (green arrows). The transition free energy 
is given by 

AG6d^6(Qi,a2) = AG/d^6(ai,a2) - AG/_6(q2). (19) 

The related transition, where a budded and dimpled domain 
coalesce to form a dimpled domain, as depicted in Fig. [U^e), 
and shown in Fig. [HIc-e) (yellow arrows) , has the transition free 
energy 

AG6d^d(ai,a2) = AG6d^6(ai, 02) - AG/^b(ai 4- 02) (20) 

+AGf^d{ai+a2). 

Comparing these two related transitions, bd ^ b and bd d, 
we see that the dominant transition is determined by the bud- 
dimple energy balance. 

The last set of transitions is when two buds coalesce to form 
a larger bud, depicted in Fig. [H^k) , with transition energy 

AG66^6(ai,a2) = AG/6-^i, (01,02) - AG/^i,(a2), (21) 

and when two buds coalesce to form a dimple, depicted in 
Fig. [9jf ) with transition energy 

AGbb^d{ai,a2) = AG/b^d(ai, 02) - AG/^6(a2). (22) 
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Comparing these two related transitions, hh ^ h and hb —> d, 
we see that the dominant transition is determined by the bud- 
dimple energy balance. 

Given the importance of the bud-dimple energy balance in 
determining the morphology resulting from a coalescence event, 
we note that if the resulting domain area is outside the range 
specified by eqn. 1101 but still larger than the critical size re- 
quired for dimpling, the dimpled morphology dominates be- 
cause the free energy change of budding is positive outside that 
range. Within this size range, selecting the dominant behav- 
ior is more subtle, and depends on the resultant domain size, 
material properties and tension. For this reason, until experi- 
mental methods are devised that can track the detailed three 
dimensional morphology of a phase separated vesicle (i.e. the 
positions and sizes of all domains and the membrane tension), 
the set of transition rules discussed in this section will remain 
largely an interpretive tool, useful for understanding the set 
of possible transitions and resultant morphologies, as well as 
their underlying physics, but difficult to quantitatively apply 
to experiment. 

We speculate that the kinetics of these coalescence transi- 
tions are either relatively fast, when diffusion is the limiting 
time scale, as might be the case in the transitions shown in 
Fig. |9ja-c,e,g,h,j,^), or relatively slow, limited by elastic inter- 
actions (Fig. [9l[d,i)) or steric hindrance (Fig. |9]^f,k)). From 
the viewpoint of coarsening of a two-phase fluid, these transi- 
tions represent new coarsening mechanisms that are linked to 
morphology, and likely have profound effects on the kinetics of 
phase separation, as demonstrated by the fact that coalescence 
of dimpled domains is inhibited by an energetic barrier. Ad- 
ditionally, these transitions suggest interesting biological pos- 
sibilities. For instance, a small volume can be encapsulated at 
a particular location, as a dimple transitions to a bud. The 
enclosed volume can then diffuse to other regions of the mem- 
brane, and either engulf more volume (see Fig. |9ji)) or deposit 
its contents at the site of another domain (see Fig. El^j)). In 
fact, both of these scenarios play out in Fig. |5][c-e). Further- 
more, it is possible that careful control of membrane tension 
[67j could regulate how large a volume is enclosed, and to which 
other domains a bud will coalesce and deposit its contents. This 
has implications for the size-selectivity of endo- and exo-cytosis 
where membrane invagination and fusion occur, as well as reg- 
ulation of plasma membrane tension |67) . 

Discussion 

The transition free energies calculated in the previous section 
have the intuitively pleasing feature of being a sum of three 
simple basis transitions (/ — > d, / 6 and ff-^f)- How- 
ever, this type of analysis is limited by the fact that it only 
admits flat, dimpled and budded as valid morphologies. More 
general theories and computational models can (and have been) 
constructed that attempt to describe all possible shapes of a do- 
main from precisely flat to fully budded, and other more com- 
plex morphologies gU [H [53l [55] [M ESI EO] . Our level of ex- 
perimental sophistication is commensurate with the simplicity 
of the analysis employed in the previous section. Conceptually, 
our model simplifles analysis by reducing domain morphology 
to one of three classes of shapes, at the cost of excluding other 
possible morphologies. Though overall an experimental minor- 



ity, domain-induced tube formation was the most common of 
these more exotic morphologies. Normally, thin lipid tubes are 
drawn out by external force [711 1721 [73] . However, in a few in- 
stances we observed domains that spontaneously collapse and 
nucleate a tube that rapidly grows many times longer than its 
persistence length, as demonstrated in Fig. 1101 Oddly, the nu- 
cleating domain is of one lipid phase, but the tube continues to 
grow from the other, majority phase by a currently unknown 
mechanism. 

In addition to limiting the class of possible morphologies, our 
analysis of morphological transitions also employs the simpli- 
fication that membrane tension is constant during a morpho- 
logical transition. In reality, our experiments take place on a 
spherical topology with constrained volume and surface area, 
such that this approximation has a range of validity. If the 
membrane area required to complete a morphological transi- 
tion is small compared to the total vesicle area (see [42| for de- 
tails) the change in membrane tension will be small. However, 
morphological transitions that require relatively large areas can 
result in significant changes to membrane tension, invalidating 
the constant tension approximation. Although, at times this 
can be an advantageous feature of our experimental system, for 
instance, when fairly small changes in vesicle area (on the or- 
der of 1%) can reduce the tension enough to cause spontaneous 
budding, as we showed earlier in this work. 

In addition to the limitations mentioned above, our exper- 
iments have a number of subtle complications. Notably, the 
task of measuring the motion and size of lipid domains is com- 
plicated by the fact that the spherical curvature of the vesicle 
slightly distorts measurements of distance and size. Addition- 
ally, the motion of domains is confined to lie in a circle defined 
by a combination of vesicle size and depth of field of the mi- 
croscope objective. We developed schemes to correct for these 
issues, as discussed in detail in the supplementary information 
of our previous work [32] • 

Summary 

Using a model multi-component membrane, we explored how 
the interplay between composition and morphology leads to 
elastic forces that spatially organize domains and significantly 
impact coalescence kinetics. We expanded upon mechanical 
models that incorporate bending stiffness, membrane tension, 
phase boundary line tension, and domain size to show that 
domains can adopt (at least) three distinct morphologies: flat, 
dimpled and budded. We showed that dimpled domains exhibit 
measurable translational and orientational order as a function 
of increasing domain areal density [491 1501 151] . Using a spheri- 
cal budding model, we showed that the transition to a budded 
state is a domain size selective process, from which one can 
estimate the membrane tension, line tension, and elastic decay 
length of a phase separated membrane. Additionally, we found 
that the large energy scales associated with changes in domain 
morphology allow us to define morphological transition rules, 
where domain size and membrane tension are likely the key 
parameters that regulate the morphological transitions. 

In the context of our understanding of the physics of phase 
separation the elastic forces between dimpled domains that 
arrest coalescence, and the morphological transitions between 
fiat, dimpled and budded domains, constitute new mechanisms 
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that govern spatial organization of domains and the temporal 
evolution of domain sizes. For cellular membranes, we specu- 
late that the elastic forces and morphological transitions can be 
controlled via careful regulation of membrane tension [67], and 
our work suggests intriguing possibilities for how small volumes 
can be encapsulated, moved, and released in a phase-separated 
membrane. 

Materials and Methods 

Giant unilamellar vesicles (GUVs) were prepared from a 
mixture of DOPC (l,2-Dioleoyl-sn-Glycero-3-Phosphocholine), 
DPPC (l,2-Dipalmitoyl-sn-Glycero-3-Phosphocholine) and 
cholesterol (Avanti Polar Lipids) (25:55;20/molar) that 
exhibits liquid- liquid phase coexistence (53]. Fluorescence 
contrast between the two lipid phases is provided by the 
rhodamine head-group labeled lipids; DOPE (1,2-Dioleoyl- 
sn-Glycero-3-Phosphoethanolamine-N- (Lissamine Rhodamine 
B Sulfonyl)) or DPPE (l,2-Dipalmitoyl-sn-Glycero-3- 
Phosphoethanolamine-N- (Lissamine Rhodamine B Sulfonyl)), 
at a molar fraction of ~ 0.005. The leaflet compositions are 
presumed symmetric and hence Co — 0. 

GUVs were formed via electroformation |23l I74j . Briefly, 
3 — 4 /ig of lipid in chloroform were deposited on an indium-tin 
oxide coated slide and dessicated for ~ 2 hrs to remove excess 
solvent. The film was then hydrated with a 100 mM sucrose 
solution and heated to ~ 50 C to be above the miscibility tran- 
sition temperature. An alternating electric field was applied; 
10 Hz for 120 minutes, 2 Hz for 50 minutes, at ~ 500Volts/m 
over ~ 2 mm. Low membrane tensions were initially achieved 
by careful osmolar balancing with sucrose (~ 100 mM) inside 
the vesicles, and glucose (~ 100 — 108 mM) outside. Using a 
custom built temperature control stage, the in situ membrane 
tension was coarsely controlled by adjusting the temperature a 
few degrees [SH [SO] . 

Domains were induced by a temperature quench and imaged 
using standard TRITC epi-fluorescence microscopy at 80x mag- 
nification with a cooled (-30 C) CCD camera (Roper Scientiflc, 
6.7 X 6.7 fim^ per pixel, 20 MHz digitization). Images were 
taken from the top or bottom of a GUV where the surface 
metric is approximately flat. Data sets contained ~ 500— 1500 
frames collected at 10-20 Hz with a varying number of do- 
mains (usually > 10). The frame rate was chosen to minimize 
exposure-time blurring of the domains, while allowing suffi- 
ciently large diffusive domain motion. Software was written to 
track the position of each well-resolved domain and calculate 
the radial distribution function. The raw radial distribution 
function was corrected for the fictitious confining potential of 
the circular geometry. The negative natural logarithm of the 
radial distribution function is the potential of mean force plus 
a constant, as shown in Fig. [3] Detailed explanations of these 
concepts can be found in the supplementary information for 

Morphological transitions were induced by quenching 
homogeneous vesicles below the de-mixing temperature and 
observing those that had many micron-sized domains. With- 
out precise control of membrane tension or the exact initial 
conditions {i.e. the exact number and size distribution of 
domains) many vesicles had to be sampled to see transitions. 
Often, a slight increase in temperature (~ 2C) was used to 



increase the available membrane area, and hence decrease the 
membrane tension enough to induce morphological transitions. 
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No. ACI-0204932, NIRT award No. CMS-0404031 and the National 
Institutes of Health award No. ROl GM084211 and the Director's 
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Figure 2: Domain morphology and coalescence, a) A nearly fully 
phase-separated vesicle, showing domains (red) flat with respect to 
the background curvature of the vesicle (blue), b) At low tension, 
domains (red) dimple and establish a non-zero boundary slope with 
respect to the curvature of the vesicle (blue), c) Flat domains on 
the surface of a vesicle - coalescence is uninhibited by elastic inter- 
actions, d) Dimpled domains on the surface of vesicle - coalescence 
is inhibited by clastic interactions between the domains, and the 
domain-size distribution is stable. Directly measuring membrane 
tension disturbs the domain size evolution, however the magnitude 
of membrane fluctuations 1471 1751 indicates that the tension in (c) is 
higher than the tension in (d). Scale bars are 10 fim. 
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Figure 3: Interactions of lipid domains as areal density increases. Left: snapshots of dimples on the surface of GU Vs. Right: corresponding 
potentials of mean force, a) At low areal density, interactions are almost purely repulsive, and there is no translational or orientational 
order — the domains are in a state analogous to a gas of particles, b-f) At higher areal density, domains 'condense' into a state where each 
domain is repelled by its neighboring domains, giving rise to energy wells that define a lattice constant and hence translational order. The 
decay envelope of these 'ringing' potentials indicates the length-scale over which the motion of domains is correlated. In all plots, the blue 
line indicates the fit to eqn. [T] where A2 is the order-correlation length and A3 is the effective lattice constant. The dashed vertical lines 
are the approximate minimum center-to-center distance between domains as determined by domain size measurement (a) or one half the 
lattice constant (b-f). Insets: Time-averaged Fourier transforms, showing that mutually repulsive clastic interactions lead to (thermally 
smeared) hexagonal order, except in (a) where the density is too low to order the domains. Scale bars are 10 /im. 
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Figure 4: Increase in correlated motion of dimpled domains as a 
function of area occupancy. This plot shows the ratio of the order- 
correlation length (A2) over the lattice constant (A3) for the vesicles 
m Fig. [3] (and one additional vesicle) as a function of the total area 
taken up by the domains divided by the total measurable vesicle 
area. The ratio A2 / A3 quantifies how many nearest neighbor domains 
{i.e. 1st, 2nd, etc.) exhibit strongly correlated motion. 




Figure 5: Shapes and energies of domain budding, a) A schematic 
of domain shape going from a fiat domain, with area a and flat 
radius p, through a dome shape with wrapping angle 8 and radius 
of curvature R, to a fully budded state, with an applied tension t. 
b) The free energy of a budding domain as a function of line tension 
(x) and wrapping angle (6) for domain size a = 10. At low line 
tension (before the blue line), both fiat and budded morphologies 
arc stable, but the flat state has a lower elastic free energy and 
there is an energy barrier between the two stable states. At the 
blue line, the free energy difference between flat and budded is zero. 
Between the blue and red lines, both morphologies are stable, but the 
budded state has a lower elastic free energy. Finally, for line tensions 
above the red line, the energy barrier disappears and budding is a 
spontaneous process. 
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Figure 6: Equilibrium phase diagram for domain budding as a 
function of dimensionless domain area and line tension. In region 
'a' flat and budded domains coexist, with flat domains at lower free 
energy. In region 'b' flat and budded domains coexist, with budded 
domains at lower free energy. In region 'c' only a single, budded 
phase is stable. The line separating regions 'a' and 'b' is given by 
Xbud/2 and between regions 'b' and 'c' by Xbud (eqn. [7]l. Dashed 
lines are trajectories of increasing membrane tension (as indicated 
by the arrows) at constant domain area. In all four trajectories 
7 = O.SfeflT/nm, Kf, = 25 kgT and tension is varied from t = 
10^^ — 10"-^ fesT/nm'^; the domain areas are Ai = •7r(100nm)^, A2 = 
7r(250nm)2, A3 = 7r(500nm)2, and A4 = 7r(1000nm)2. 
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Figure 7: Size-selective domain budding, a) Dimples on the sur- 
face of a GUV are initially arranged by their repulsive interactions, 
b) and c) A slight increase in temperature decreases membrane 
tension, causing smaller domains to spontaneously bud (marked 
with red arrows) and wander freely on the vesicle surface, while 
larger domains remain dimpled. The mean size of budding domains 
is r ~ 0.93 ± 0.18 /^m from which we estimate a line tension of 
7 ~ 0.45 kgT/mn. d) Plot of the potential of mean force between the 
dimpled domains in (a), moments before inducing spontaneous bud- 
ding, e) Budding stability diagram, showing solutions to dG/d8 = 
with 7 = 0.45 kgT/raa. Solid red lines are stable solutions at energy 
minima; dashed red lines are unstable solutions on the energy bar- 
riers. Regions 1 and 3 are coexistence regimes, while region 2 is a 
spontaneous budding regime, only stable at \0\ = it. The blue dots 
indicate domain areas with radii r ~ 0.93, 0.93 it 0.18 fira. f) The red 
curve shows the free energy of budding at r = 1.2 X 10"'^ fc^T/nm^, 
which is greater than zero for all domain sizes, and hence all do- 
mains would remain flat/dimpled. The green curve shows the free 
energy of budding at r = 2.4 X 10"'* kgT/nm^ . Domains with radius 
r ~ 0.25 — 3.5 iim have a negative free energy of budding, all other 
sizes remain flat/dimpled. Most domains within this size range must 
still overcome an energy barrier to bud, but for a small range of do- 
main sizes (r ~ 0.75 — 1.11 lira), indicated by the blue line segment, 
budding is a spontaneous process. The energies are calculated using 
Kb = 25 kgT and 7 ~ 0.45fcBr/nm. In (a-c) scale bars are 10 /im. 
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Figure 8: Gallery of morphological transitions, a) Two dimpled domains (indicated by the white arrows) interact on the surface of vesicle, 
eventually coalescing to yield a larger dimpled domain (sec Fig. |9ld)). b) An equatorial view of a dimple-to-bud transition (indicated by 
the red arrows), c-e) Time courses of multiple types of morphological transitions. Arrows are color-coded and point to before and after 
each transition: red arrows indicate a dimple to bud transition, green arrows indicate a bud engulfing a dimple to form a larger bud (see 
Fig. Hj)), and yello w arrows indicate a bud recombining with a larger dimpled domain (see Fig. |9le)). Using video microscopy, we can 
put an upper bound on the time scale of the d — > 6, dft — > d and dft — > 6 transitions at ~ 200 ± 80 ms, ~ 160 it 70 ms and ~ 210 it 70 ms, 
respectively. Scale bars are 10 ^m. 
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Figure 9: Algebra of Morphology. All rows show additive free energies of transition from an initial state on the left to a final state on 
the right, a) Two flat domains coalesce via diffusion, yielding a flat domain, b) Domains, too small to dimple, coalesce to attain a size 
capable of dimpling, c) A dimple and a domain too small to dimple coalesce to yield a larger dimpled domain, d) Two interacting dimples 
coalescence, yielding a larger dimpled domain, e) A dimple and a bud coalesce to yield a larger dimpled domain, f) Two buds coalesce to 
yield a dimpled domain, g) A flat domain coalesces with a bud, yielding a dimpled domain, h) A flat domain coalesces with a dimpled 
domain to yield a bud. i) Two interacting dimples coalescence, yielding a bud. j) A bud coalesces with a dimple, yielding a larger budded 
domain, k) Two budded domains coalesce to form a larger budded domain. 1) A flat domain coalesces with a bud to yield a larger budded 
domain. 
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Figure 10: Domain-nucleated spontaneous tube formation. A time 
series of spontaneous tube formation, nucleated from a domain (as 
indicated by the white arrow). This relatively uncommon morphol- 
ogy is not explained within the context of our simple model. The 
lipid tube (bright) is many times longer than its persistence length, 
yet perplcxingly, grows from the tube tip. With limited optical res- 
olution, we estimate the tube diameter to be < 500 nm. The scale 
bar is 10 /im. 
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